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An octupole ordering model is studied by the mean field theory, and its relevance to the 
phase IV of Ce^Lai-^Be is discussed. The observed lattice distortion along the [Iff] direction 
is interpreted in terms of the Tsg-type /erro-quadrupole moment induced by an antiferro- 
octupole ordered state with r$ u symmetry. The octupole model also accounts for the cusp in 
the magnetization as in the Neel transition, and the softening of the elastic constant C44 below 
the ordering temperature. However, the internal magnetic field due to the octupole moment is 
smaller than the observed one by an order of magnitude. Also discussed is the possibility of a 
pressure induced antiferromagnetic moment in the octupole-ordered state. 

KEYWORDS: rare-earth hexaboride, Ce^Lai-^Be, NdE>6, NpC>2, octupole ordering, elastic constant, ther- 
mal expansion, internal field, pressure induced moment 
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1. Introduction 

Orbital orderings in /-electron systems have attracted 
much attention. In usual cases, an orbital-ordered phase 
is a quadrupole-ordered phase. However, higher multi- 
pole orderings can also be realized. Indeed, possibilities 
of octupole orderings in Np02 1 ~ 4 - ) and in Ce^Lai-^Bg 5, 6 ) 
have been discussed. In this paper, we study an octupole 
ordering model for the phase IV of Ce^Lai-^Bg paying 
attention to the coupling between /-electrons and lattice. 

CeB 6 has a cubic crystal structure of the CaB 6 -type, 
and Ce ions form a simple cubic lattice. In this ma- 
terial, there are three phases; the paramagnetic phase 
(called phase I), the antiferro-quadrupole phase (phase 
II) and the antiferromagnetic phase (phase III). The role 
of the octupole interaction in these phases has been dis- 
cussed 6-9 - 1 to explain experimental observations. In the 
dilute alloys Ce^Lai-^Bg, the so-called phase IV was 
found at x ~ 0.75. 10 -* This phase has attracted much at- 
tention. A large softening of the elastic constant C44 just 
below the transition temperature Tt_iv from the phase I 
to the phase IV was observed. 11 '' Furthermore, the mag- 
netic susceptibility shows a cusp at T\.\\, and the mag- 
netization is almost isotropic at ambient pressure. 12 -* 

It was recently revealed that the lattice shrinks along 
the [111] direction in the phase IV. 13 ' 14 ) It is probable 
that the softening of C44 is due to this distortion. It 
may be tempting to ascribe the distortion to the ferro 
r5g-type quadrupole order. However, the quadrupolc mo- 
ment is not necessarily the primary order parameter. 
Indeed, the temperature dependence of C44 is different 
from usual ferro-quadrupole-ordered materials. 11 ) In the 
phase I of Ceo.75Lao.25B6, the temperature dependence 
of C44 follows the Curie- Weiss law with an antiferro- 
quadrupole interaction, but softens drastically just below 
21-iv- The entropy changes more at the transition from 
the phase I to the phase IV than at the other transition 
from the phase IV to the phase III. 11 ' 15 ) On the contrary, 
the entropy changes very little at the quadrupole tran- 
sition in CeBg. 15 -* Thus it is likely that the degeneracy 
of each /-electron state is already lifted in the phase IV. 



This fact also makes it difficult for the T 5g quadrupoles 
to be the primary order parameter in the phase IV. 

A large change in the internal field at Tt-ivj as probed 
by NMR 16 ) and /iSR, 17 ' suggests strongly that the time 
reversal symmetry is broken in the phase IV. A pure 
quadrupole order is incompatible with the broken time 
reversal symmetry. In addition, a neutron diffraction ex- 
periment found no magnetic reflection in the phase IV. 18 ) 
Thus, dipole moments are also unlikely to be the primary 
order parameter, although the time reversal symmetry is 
broken. Therefore, octupole moments, which break the 
time reversal symmetry, become a candidate for the order 
parameter in the phase IV. 5, 6 ) In addition, the order pa- 
rameter should have an anisotropic nature, because the 
cubic symmetry is broken, and the anisotropy in the mag- 
netization develops under uniaxial pressure. 19 ~ 21 ) Thus 
the r 5tl -type octupole moment, among all octupole mo- 
ments, is the most plausible candidate for the order pa- 
rameter in the phase IV. We will discuss this point in 
detail in § 2. 

Kusunose and Kuramoto 6 ^ have already pointed out 
using the Ginzburg-Landau (GL) theory that the T^ u 
octupole order with a finite wave number should accom- 
pany a ferro-quadrupole moment, and have suggested 
a possible lattice distortion. However, evaluation of the 
magnitude of the distortion is beyond their GL theory, 
because the GL theory is appropriate near the transition 
temperature. In this paper, we explore in much greater 
detail the consequence of the T 5u octupole order by the 
mean field theory, and propose that the lattice distortion 
and the softening of the elastic constant C44 are due to 
the order of the T§ u octupole moment. Some of the re- 
sults in this paper were reported in a recent paper. 22 ' In 
this paper, we describe the details of the calculation, and 
add some new results, especially, the quadrupole suscep- 
tibilities. 

This paper is organized as follows. In § 2, we introduce 
an octupole ordering model and solve this model by using 
the mean field theory. In § 3, § 4 and § 5, we discuss the 
magnetization, elastic constants and lattice distortion, 
respectively, in the octupole-ordered state. In § 6, we 
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estimate the internal field from the octupole moment. In 
§ 7, we discuss a possible dipole moment under uniaxial 
pressure in the octupole-ordered state. The last section 
is devoted to a summary. 

2. Mean Field Theory 

The crystalline electric field (CEF) ground state of 
Ce 3+ (J = 5/2) in Ce a; Lai_ x B 6 is the T 8 quartet. 23 - 25 ) 
The excited level 1^ lies about 500 K from the r 8 level 
and is neglected. The r 8 states are represented in terms 
of cigenstates of J z as 



T> 

I) 
T> 




'5/6| - 5/2) 
1/2), 



3/2), 
3/2), 



I- I) = 1-1/2). 



(2.1) 
(2.2) 
(2-3) 
(2.4) 



Within the r 8 quartet, the number of independent multi- 
polar moments is 15, and active octupole moments have 
cither T2 U , ^4u or T^ u symmetry. 7 ) 

We first discuss the properties of the T 2u and r 4u 
moments. The r 4 „-type octupole moments accompany 
dipole moments, because dipole moments have the same 
symmetry T^. 7 ) Thus the T^ u octupole moments are 
unlikely to be the order parameter in the phase IV. The 
r2«-type octupole moment, as proposed in ref. 5, is pseu- 
doscalar, and eigenstates of the T 2 u moment have the cu- 
bic symmetry. Hence it does not accompany quadrupole 
moments, and it seems difficult to explain the lattice dis- 
tortion in the phase IV by the T 2u octupole orderings. In 
the followings, we study the properties of r 5 „ octupole 
moments, and discuss their possibility as the order pa- 
rameter of the phase IV. 

Let us introduce pseudospins <t and r to describe the 
T 8 quartet: 



T z \±a z ) =±|±ct z ), 



(2.5) 



t 2 |t 2 T)=+|t*T), ° z \t z l) = -\r z |), (2.6) 

and the transverse components which flip the pseu- 
dospins. The relevant multipolc moments are given with 
the notation (a, (3, 7) = (x, y, z), (y, z, x) or (z, x, y) by 
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and 



(2.7) 
(2.8) 
(2-9) 
(2.10) 
(2.H) 

(2-12) 



where bars on the products represent normalized sym- 
metrization, e.g., J X J 2 = (J x Jy + J y J x Jy + JyJ x )/i, and 



A 



5g- 



r 8 



'V3' 
r 5g quadrupole 
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(«V<0) 
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-4lA iu 



T 5u octupole 



Fig. 1. The level scheme in the T$ u interaction (right) and in 
the interaction (left). The quantities A 5u and A 5g are the 
octupolar and quadrupolar mean fields. 



we have introduced the notations 

C ± = -\(t x ±V3t z ). 



(2.13) 
(2-14) 



2.1 Level Scheme 

Before introducing a octupole ordering model on 
a lattice, we discuss the 'easy axis' of the r 5 „ 
octupole moment under a fictitious octupole field 
A 5u = (Al u ,Al u ,A 5 z u ). It is sufficient to consider high- 
symmetry directions of A bu . With |A 5 "| = 1, the four 
eigenvalues (Ev) of relevant cases are given as follows: 

(i) A 5u || (0,0,1): Ev[T 5 "] = ±1 (both doubly degener- 
ate); 

(ii) A 5u || (1,1,0): Ev[(T^ + r 5 «)/V2] =±(>/3±l)/2; 

(iii) A 5u || (1, 1, 1) : Ev[(T 5u + + T z 5 «)/V3] = ±V2 
and (doubly degenerate). 

Case (iii) gives the largest eigenvalue, and thus the easy 
axis is along the [111] direction and equivalent ones. 

In the mean field theory, the fictitious octupole 
field represents the intersite interaction. When 
the intersite interaction is isotropic, a state with 

(\( T x U }\A( T y U )\A( T ! u )\) II (1,1,1) is the most stable. 
We note that the operator T 5tl + T^ u + T 5tl commutes 
with O yz + O zx + O xy , and these two operators can 
be diagonalizcd simultaneously. On the right side of 
Fig. 1, we illustrate the level scheme in octupole field 
A 5u = A 5u (l,l,l)/V3. The T 8 level splits into three 
levels. Not only the time reversal symmetry, but also 
the cubic symmetry is broken. The multipolar moments 
in the ground and highest states are 



(T x 5 "> = (T 5 «) - (T 5 «) = tVV 
(O yz ) = (O zx ) = (O xy ) = -2/3, 



(2-15) 
(2-16) 



and the others are zero. Thus, ferro, antiferro and other 
collinearly ordered states with ((T 5tI ), (T 5u ), (T 5tt )) || 
(1,1,1) have a homogeneous T^g moment. With the 
quadrupole-strain coupling, the crystal should distort 
along [111]. 

It is obvious that a r 59 -type ferro-quadrupole inter- 
action alone can lead to (O yz ) = (O zx ) = (O xy ) 7^ 0. 
The Kramers degeneracy remains in this case. The en- 
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Fig. 2. Temperature dependence of the antifcrro-octupole mo- 
ment (T z 5u > (= (Tj? 11 ) = (Ty 11 )). The dashed line is the asymp- 
totic behavior V2(l - T/T^) 1 / 2 at T ~ T 5u . 



Fig. 3. Temperature dependence of the ferro-quadrupole moment 
(Oxy) (= (Oyz) = (O zx ))- The dashed line is the asymptotic 
behavior -2(1 - T/T 5u ) at T ~ T 5u . 



ergy level splits into two levels, and the ground state is 
twofold degenerate as shown on the left side of Fig. 1. 

2.2 Antiferro-Octupole Ordering Model 

We proceed to consider the Ts quartets on a simple 
cubic lattice. From the above discussion, we find that 
the order parameter which breaks the time reversal sym- 
metry and induces quadrupole moments, but which ac- 
companies no dipole moment is only the T^ u octupole 
moment. However, we cannot determine the periodic- 
ity of the r 5u -ordered state in the phase IV only from 
the above consideration, because any collinearly ordered 
state with ((T^),(T^),(T^)) || (1,1,1) has a uni- 
form T 5g moment. In the present study, we consider a 
G-type (staggered) antiferro-octupole order as the sim- 
plest example of T^ u orders. The importance of the T^ u 
nearest-neighbor interaction in causing the change from 
the phase III to the phase III', even with a weak mag- 
netic field, was pointed out in ref. 6. We consider only 
this nearest-neighbor interaction, and take the following 
model: 



n = i 5u J2 



E 



(2.17) 



where denotes a nearest-neighbor pair. We study 

this Hamiltonian by the mean field theory, and choose the 
value of I 5u so as to reproduce the transition temperature 
Thv in Ceo.75Lao.25Be, i.e., T 5u = 6I 5u = 1.7 K. In 
carrying out the mean field theory, we assume the G-type 
two sublattice structure, but do not assume a particular 
solution, e.g., (T^ w ) = (T% u ) = (T z 5m ). When we obtain 
several self-consistent solutions, we choose the one which 
has the lowest free energy. 

The temperature dependences of (T^ u ) and (O xy ) are 
shown in Figs. 2 and 3. The solution obtained has the 
antiferro-octupole order with (T| M ) = (T^ u ) = (T z 5u ), 
accompanying the ferro-quadrupole moment (O yz ) = 
(Ozx) = (O xy ). At T ~ T 5 „, we obtain (Xf 1 ) oc (1 - 
T/T^u) 1 / 2 as a consequence of the mean field theory. 
On the other hand, the quadrupole moment behaves as 
(0 X y) oc (1 — T/T^u), which indicates the feature of the 
induced order parameter. 6 ) We also find other equivalent 
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Fig. 4. Temperature dependence of the magnetization in mag- 
netic field H = 0.2 T along various directions. The domains are 
selected by the magnetic field directions. 



solutions: ((T Q 5 «), (T|«), (T^), (O 0J ), (O ja ), (O afj )) = 
(±B,±B, T B,+C,+C,-C) with (a,/3, 7 ) = (x,y,z), 
(y 7 z 7 x) or (z,x,y), where quantities B and C are posi- 
tive and depend on temperature. The degeneracy of these 
solutions can be lifted by a magnetic field or pressure. In 
the following sections, we choose a domain in which the 
free energy is minimized. 

3. Magnetization 

A magnetic field can lift the degeneracy of the do- 
mains. For H [111], the (111) domain given by 
cqs. (2.15) and (2.16) is not the most stable one. In- 
stead, a state where two of (0 Qj g)'s have positive val- 
ues, and the other has a negative value is stabilized. For 
H || [110], (J x ) and (J y ) have positive values, and a 
state with (O xy ) > is realized. On the other hand, all 
the domains obtained in § 2 have the same energy under 
H || [001]. 

In Fig. 4, we show the temperature dependence of 
the magnetization in magnetic field H = 0.2 T along 
the three high-symmetry directions. The magnetization 
has a cusp at the transition temperature T 5u . The pres- 
ence of a cusp is consistent with experimental observa- 
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0.5 1 1.5 2 2.5 

r(K) 

Fig. 5. Temperature dependence of the magnetization in the sin- 
gle domain with (T|») = (T^ u ) = {T% u ) in magnetic field 
H = 0.2 T along various directions. 

tions. 12, 19 ~ 21 ) Very recently, magnetization under uni- 
axial pressure p = 100 MPa along the [111] direction 
has been measured at H = 0.2 T and H = 0.4 T. 21 ) 
The observed magnetization is anisotropic and Mni i > 
Mjooi] > Afhiji. It seems that the system under the uni- 
axial pressure has a single domain. In order to deal with 
this situation, we have computed the temperature de- 
pendence of the magnetization assuming a single domain 
with (T^) = (2*«) = (2*«) ((O yz ) = (O zx ) = (O xy )) 
in magnetic field H = 0.2 T. The results are shown in 
Fig. 5. The anisotropy derived by our model is consistent 
with the experimental observation. 

On the other hand, a previous experiment without 
an external stress, the magnetization below Tj_iv is al- 
most isotropic. 12 > The observed isotropy may be due to 
a multi-domain structure. Our model have four equiv- 
alent solutions in the absence of a magnetic field. If 
these four domains are distributed with the same prob- 
ability, the magnetic susceptibility is isotropic, because 
the cubic symmetry is effectively recovered in the multi- 
domain state. The anisotropy in the magnetization under 
p || [001] was also reported. 19 ' 20 ) It is not easy to explain 
this anisotropy, because the four domains in our model 
are equivalent under p || [001] and it seems that a single 
domain state is not realized under p || [001]. 

We have also calculated the magnetization in the case 
of ferro-octupole interaction. The result is qualitatively 
similar to that in the antiferro-octupole case; the mag- 
netization has a cusp, and the anisotropy in the single 
domain with (T^ u ) = (T^ u ) = (T^ u ) is the same as that 
of the antiferro-octupole case. Therefore, it is difficult to 
exclude the possibility of the ferro-octupole ordering in 
the phase IV from the magnetization alone. 

For comparison, we show the magnetization of a model 
in which only the T^g ferro-quadrupole interaction is 
present in Fig. 6. The interaction between quadrupole 
moments is chosen so as to reproduce the transition tem- 
perature Ti-iv = 1-7 K. In this calculation, the domains 
are determined by the direction of a magnetic field. We 
cannot fix a domain in this case, because the ground state 
is continuously degenerate in the absence of a magnetic 
field. For example, the state with (O yz ) = (O zx ) = {O xy } 




T(K) 

Fig. 6. Temperature dependence of the magnetization in mag- 
netic field H = 0.2 T along various directions for the case where 
the quadrupole moment is the primary order parameter. 

is not even metastable in H || [001]. The magnetization 
only slightly changes at the quadrupole transition tem- 
perature, and does not have a cusp. Thus it is unlikely 
that the primary order parameter in the phase IV is the 
ferro-quadrupole moment. 

4. Elastic Constant 

In ferro-quadrupole-ordered materials such as 
TmZn, 26 ) the elastic constant often shows a large 
softening, following the Curie- Weiss law, as the transi- 
tion temperature is approached from above. However, 
the softening of C44 above Ti_iv is very small in 
Ceo.75Lao.25B6, and a large softening is observed just 
below Tt.p/. 11 -' Moreover, an antiferro-quadrupole in- 
teraction g' T = — 2 K is deduced by a fitting of the 
temperature dependence of C44 above Ti.py Thus it is 
difficult to ascribe this softening to a ferro-quadrupole 
interaction. In this section, we consider the elastic 
constant in the octupole ordering model. We do not 
include the quadrupolc-quadrupole interaction g' r& for 
simplicity. 

The elastic energy associated with the quadrupole mo- 
ments is given per unit volume by 

E = E (rrA 0) +SrE £ r 7 Or> j , (4.18) 
r 7 V i J 

where er 7 is the strain tensor, gr is the magneto-elastic 
coupling constant, and C r °' is the bare elastic constant. 
The second sum runs over Ce sites i in a unit volume. 
The second derivative of the free energy with respect to 
ep 7 gives the elastic constant Cr, 27 -* which is given by 

C T = 4 0) - Nglxr, (4.19) 
where N is the number of Ce ions in the unit volume and 
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Fig. 7. Temperature dependence of the quadrupolc susceptibility 
Xr 3 (solid line). The Curie law behavior \r 3 = 64/(3fc B T) at 
high temperatures is also shown (dashed line). 



Fig. 8. Temperature dependence of the quadrupolc susceptibility 
Xr 5 (solid line). The Curie law behavior Xp = 4/(3£;bT) at high 
temperatures is also shown (dashed line). 



Xt is the quadrupole susceptibility denned by 
d(Or 7 ) 

Srer T =0 

dE n 



XT 



d (.9r£r 7 ) 



1 



^2p n [{n\0 Tl \n) - (Or 7 )o] 



d(.9rer 7 ) 



d(n|Or 7 |^} 



4^ Pn "d( 5 rer 7 ) 



3r<r T =0 



9rer-,=0 



(4.20) 



where E n and \n) are the eigenvalue and eigenstate of the 
total Hamiltonian, and (■■■)() denotes the expectation 
value without the strain. The occupation probability p n 
is given by 

e -E n /k B T 

Pn = — -E m/ks T - ( 421 ) 



The elastic constants, strain tensors and quadrupole mo- 
ments with r 3 and T 5 symmetries are given as follows: 



Cr 3 



(Cu - C 12 )/2, (4.22) 
4C 44 , (4.23) 
e„ = (2e zz - e xx - e yy )/V3, (4.24) 

£y ^xx ^yyi (4.25) 

(4.26) 
(4.27) 
(4.28) 

O aP . (4.29) 

We calculate xr by differentiating (Or 7 ) with respect to 
gv^v-i numerically. In Figs. 7 and 8, we show the temper- 
ature dependence of xr 3 and xr 5 ■ Above the transition 
temperature T 5u , the octupole interaction I 5u does not 
influence xr 3 and xr 5 , because the strain er 7 does not 
induce octupole moments in the non-ordered state. The 
quadrupole susceptibility xr 5 jumps at T 5u due to the 



£ r 3 tt 
er 3 v 
e r 6 a/3 
Or 3U 
Or 3V 



e af3, 
(1° 

n 2 



o 



r 5 a/3 



octupole order. Although the mean field theory gives a 
finite jump at T bu , fluctuations beyond the mean field 
theory should lead to divergence of xr 5 at T bu . 

For further understanding of the jump obtained above, 
we explore xr 5 at T ~ T 5u analytically. The mean field 
Hamiltonian for the A-sublattice with finite e xy without 
magnetic field is given by 

Hmf = -6I 5u J2( T " U ) T t + •'/> < ° • ( 4 -30) 

a 

For T < T 5u and e xy ? 0, we have (T|») = (T^) ? 
(T z u ). We determine these quantities self-consistently. 
We obtain E n and \n) of 7^mf by perturbation theory 
and calculate xr s by using eq. (4.20). At T ~ T 5u , we 
obtain 



XT:, = 



4 2VQI 5u 
3fceT k^T 



Xx u - xy h, (4.31) 



where 



knT 



T 5« + T 5^ 0) (4 32) 



with Eg the mean field ground state energy per Ce ion. 
We have also introduced the octupole susceptibility with 
respect to the strain e xy , which is defined by 



X 



bu-xy 



d(gr 5 e xy ) 



gr 5 t 



By symmetry, we obtain 



5u-xy _ bu-xy 



(4.33) 



(4.34) 



At T > T bu , we have h = and xl u ' xy = 0, and then only 
the first term in eq. (4.31), i.e, the Curie term remains 
finite. The second term in eq. (4.31) becomes finite at 
T < T 5u . The octupole susceptibility satisfies the equa- 
tion 



h 2 



y5u-xy 

■Y^u-xy 
Ax 



fn/6 
k B T 



(4.35) 
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X 



5u-xy 



When h ^ 0, i.c, T < T$ u , we obtain 



2 AS 1 

1 4V6 
k B T 3 



X. 



5u-xy 



h-\ 



(4.36) 



(4.37) 



Therefore Xx X!/ diverges at Ti_rv, and the second term 
in eq. (4.31) becomes finite just below T$ u . The change 
of the susceptibility Axr 5 at T 5 „ is given by 



Axr 5 = 



(4.38) 



3fcBT 5 „ 

Experimentally, C44 remains finite at Ti_rv- Moreover, 
C44 becomes smaller as temperature is lowered below 
Ti-iv- In general, the second-order transition involving 
the lattice distortion along [111] leads to divergence of 
C44 at the transition temperature. Thus the observed 
behavior of the elastic constant seems to be influenced 
by the experimental setup, and may not be intrinsic 
to C44. One possible reason for this behavior is that 
the large ultrasonic absorption accompanying softening 
makes it difficult to determine the elastic constant ac- 
curately. Another possibility is the mode mixing of the 
sound wave. To determine C44, sound wave with wave 
vector k || [001] and polarization u j [100] was used 
in the experiment. 11 ' However, this mode is not a nor- 
mal mode below Ti_rv, because the cubic symmetry is 
broken. It is desirable to measure the longitudinal ul- 
trasound velocity with k || [111] which is expected to 
reflect the anomaly in C44 directly, because this mode 
remains a normal mode even in the trigonal symmetry 
((O yz ) = (O zx ) = (O xy ) ? 0). 

5. Lattice Distortion 

We now discuss the lattice distortion due to the 
antifcrro-octupolc order. By minimizing the elastic en- 
ergy (4.18) associated with the T^ g moments, we obtain 

9r 5 



4C 



(5.39) 



We use the following experimental values for 
Ceo.75Lao.25Be: |flr B | = 155 K, CH } ~ 8.2 x 
10 11 erg/cm 3 , 11 ^ and the lattice constant a = 4.13 A. 
At absolute zero, the magnitude of quadrupole mo- 
ments, induced by the octupole order, is given by 
\(O yz i) \ = \{O zx i) \ = \{O xy i) \ = 2/3, and we obtain 

|e„*| = = kxvl = 4.6 x 1(T B . (5.40) 

In order to account for the observed lattice contraction 
along [111], 13 ' 14 ) we consider two possibilities: (i) gr s < 
and (ii) gr s > 0. 

5.1 case (i) <?p 5 < 

In case (i), we assume that a small stress which accom- 
panies the measurement breaks the equivalence of four 
octupole domains, and chooses a domain for which the 
contraction becomes maximum along [111] with <?r s < 0. 
In this case, the domain with (O yz i) — (O zx ,) = 
(O xy i) = —2/3 is chosen, and all e a p are negative. Then 



we obtain 

12 1 
AZ/Z = -e B + - (e yz + e zx + e xy ) = -e B + 2e xy , (5.41) 

along [111], where e B is the volume strain which is not 
included in our model. Experimentally, the value of the 
volume strain in Ceo.7Lao.3B6 was obtained to be 6b = 
8x 10~ 6 at 1.3 K. 14 ) Note that along [111], [111] and [111] 
directions, the lattice should expand by 



1 2 

AZ/Z = -e B + - 1 esj/l 



(5.42) 



5.2 case (ii) gr 5 > 

In case (ii), a positive stress along [111] may favor a 
domain which contracts along this direction. In this case, 
two of e Q/ 3 are negative and the other is positive. Then 
the contraction along [111] is given by, 



1 



1 



AZ/Z = -e B + -={e yz + e zx + e xv ) = -e B 



g-ij ■ g \-y* 1 ■ ~u,y/ g~±J (5-43) 

The crystal expands along another direction. In the single 
domain with (O yz ) = (O zx ) = —(O xy ), for example, an 
expansion along [111] should be present, and it is given 

by 

12 1 

AZ/Z = -CB + ^(-£ yz -e zx + e xy ) = -e B + 2\e xy \. (5.44) 

5.3 Comparison with the Experiment 

According to an experiment, 14 ) the shear strain on the 
assumption of the trigonal symmetry around [111], i.e., 
case (i), is derived as e a p = — 4x 10~ 6 without a magnetic 
field at 1.3 K. Assuming that the phase IV is stable down 
to zero temperature, e a p is extrapolated to —6 x 10~ 6 
to —10 x 10 -6 at absolute zero. The absolute value is 
by an order of magnitude smaller than our estimate in 
cq. (5.40). Thus the observed reduction of |(O a| g)| should 
come from quantum fluctuations neglected here. 

Another possibility is that case (ii) is realized. In this 
case, the value of the strain becomes three times larger, 
i.e., \e a fj\ = (2-3) x 10~ 5 , when deduced from the same 
experimental result. 14 ^ This value is comparable to that 
given in cq. (5.40). As mentioned in § 3, a state where 
one of (O a /3)'s has a sign different from the others is 
stabilized under H \\ [111] without an external stress. 
Experimentally the contraction is enhanced under H | 
[111]. 13 ' 14 ) This feature favors case (ii) in our model. 

To determine which case is realized in the experiment, 
it is desirable to measure the lattice distortion along the 
[110] direction. The lattice distortion along [110] is given 
by AZ/Z = (e B /3) + e xy . Thus, in case (i), our model pre- 
dicts a smaller magnitude of lattice contraction along the 
[110] direction as compared with the contraction (5.41) 
along the [111] direction. On the other hand, in case (ii), 
the magnitude of lattice contraction along [110] is ex- 
pected to be larger than that of the contraction (5.43) 
along [111]. 

54 NdB 6 

For comparison, we estimate in the same manner the 
magnitude of the lattice distortion of NdB 6 . This ma- 
terial undergoes an A-type antifcrromagnctic order at 
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Tn — 8 K 28 ) accompanying a lattice distortion along 
the [001] direction. 29 ) This lattice distortion indicates 
a 0° ferro-quadrupole order. We estimate the lattice 
distortion in the OS, ferro-quadrupole-ordered state by 
using the experimental values: \gr 3 \ — 220 K, (Cn — 
Ci 2 )/2 ~ 20 x 10 11 erg/cm 3 , 30 ) a = 4.12 A, and the 
Lea-Leask- Wolf parameter xllw = — 0.82. 31 ) We assume 
that (C^ - C$)/2 is almost the same as the observed 
(Cn — C\i)j2. In a magnetic field H || [100], a state with 
the antifcrromagnetic moment perpendicular to [100] is 
stabilized. For a domain where the magnetic moment is 
along [001] or [010], we obtain 



\Al\/l = 2.8 x 10" 



(5.45) 



along H, on the assumption e B = 0. This value is fa- 
vorably comparable to the experimental one Al/l ~ 
-2 x 10~ 4 at 2 K and H = 2.1 T. 29 ) The distortion 
in NdB 6 is by an order of magnitude larger than those in 
eqs. (5.41) and (5.43). The reason is that the quadrupole 
moment 1(0°) I = 4-5 in NdB@ is much larger than the 
corresponding value \ (O xy }\ = 2/3 in the octupolar state 
of Ce^Lai-zBg, while the magneto-elastic coupling con- 
stants are of the same order. 

6. Internal Field 

In this section, we discuss the internal magnetic field 
associated with the octupole order. The fiSR time spec- 
tra in the phase IV consist of a Gaussian component 
and an exponential component. 17 -* The observation of 
a Gaussian relaxation indicates that internal fields are 
randomly distributed and/or fluctuating. The internal 
field deduced from the Gaussian relaxation is of the or- 
der of 0.1 T. We discuss whether the octupole moment 
can be the origin of the Gaussian relaxation. In the /iSK 
measurement, [i + locates at (a/2,0,0) and equivalent 
sites 32-34 ' with a Ce ion chosen as the origin. As a ref- 
erence the internal field from a Bohr magneton [is is 
estimated to be 

#di P oic = ^2)3 <" B — °' 1T ' (6.46) 

with a/2 ~ 2 A. The internal field from an octupole with 
the size r is estimated to be 

J2. 



octupole 



(a/2) 



-Mb ^ 0.01T, 



(6.47) 



with r ~ as — 0.53 A. For a more accurate estimate, we 
consider the multipole expansion of the vector potential 
from local electrons, which is given by 35 -* 

A ^ = E T r ~ (fe+1) ( zc, ™ )(0 ' ^) Mr ' (6 ' 48) 

k,m 

where I is the orbital angular momentum operator, 
Cm' (6, <ft) is y/ 47r/ (2k + 1) times the spherical harmonics 
Yk m (Q,4>), and M™ is the magnetic multipole moment. 
The multipole moment is determined by the wave func- 
tion ipi(r) of the i-th electron, and the orbital and spin 



angular momentum operators /, and s». Namely we have 

Mr= M £/d r .*nn>(v4- c <«-(».,*)) 



(6.49) 



Equation (6.49) is evaluated with the use of the operator 
equivalents method. 36 ) For our purpose, it is sufficient to 
consider only octupole moments. Then we obtain M™ 
through calculation of the reduced matrix element of the 
third-rank tensor. The result for one-electron states with 
J = 5/2, L = 3, S = 1/2 is given by 

2 



MS 



35 



Mr 2 >(Ji 3) >, 



(6.50) 



where the third-rank tensor operators jffl are defined by 



^3 ' — . (Jx + iJy) 3 , 



J. 



(3) 



1 



y/(3 + m)(4 - m) 
We introduce linear combinations: 



(6.51) 

[J x -U y ,J^]. (6.52) 



(6.53) 



Mt x = \l-V5(M+ 3 - M 3 - 3 ) + V3(M 3 +1 - M3- 1 )], 

(6.54) 

Mt" = ^[V5(M 3 +3 + M 3 - 3 ) + V3(M 3 + 1 + M3- 1 )], 

(6.55) 



Ml uz = M 3 °, 



(6.56) 



5«x = I[V3(M+ 3 - M 3 - 3 ) + V5(M 3 +1 - M^ 1 )], 



M, 

(6.57) 

Ml uy = ^[V3(M+ 3 + M 3 - 3 ) - ^(M 3 +1 + M 3 '% 

(6.58) 

M |- = i=[M 3 +2 +M 3 - 2 ]. (6.59) 
We then obtain H = V x A as 
H z (r) = J-^ [6Vl50 xy (r) (7V30 2 °(r) + 4) M, 

+150 zx (r) (7V30 2 te (r)-2)M ; 
+\W yz (r) (7V30° y (r) - 2) M^ y 
+2 {\{RO\(r)0%(r) - 20\/3O 2 ) (r) - 28) M- 
+9Vl50 zx (r) (70 2 » + 2) 
+9VTbO yz (r) (j0 2 2 y (r) - 2) M^ v 

+6VT50^(r) (jV30%(r) + 4) M 3 5 " z ] , 

(6.60) 



j4ux 



4u Z 
3 



8 J. Phys. Soc. Jpn. 



Full Paper 



Katsunori Kubo and Yoshio KuRAMOTO 



where 

Of{r) = {2x 2 - y 2 - z 2 )/(V3r 2 ), (6.61) 

O 2 y (r) = (2y 2 -z 2 -x 2 )/(V3r 2 ), (6.62) 

2 2 x (r) = (y 2 -z 2 )/r 2 , (6.63) 

2 2 y (r) = (z 2 -x 2 )/r 2 , (6.64) 

and the other 0(r)'s are obtained by replacing J a in 
eqs. (2.8)-(2.10) with a/r. The other components H x (r) 
and Hy(r) are obtained by exchanging x,y,z cyclically 
from above. In the Tsu-ordered state, we obtain M^ ux — 
Ml uv = Ml uz ^ 0, and eqs. (6.53)-(6.56) are zero. 

To estimate the internal magnetic field, we use the 
value (r 2 ) = 1.298 in atomic unit, which was obtained 
with a relativistic Dirac-Fock calculation by Freeman and 
Desclaux. 37 ) Then we obtain about 40 G as the internal 
field at (a/2,0,0) from a Ce ion. This value is by an 
order of magnitude smaller than that derived by the /iSR 
measurement. Thus the static octupole moment alone 
cannot account for the Gaussian relaxation of the fiSR 
spectra. In a future work, we plan to study in more detail 
fluctuations in the octupole-ordered state. If the origin 
of the fj,SK spectra is not an effect of fluctuations, we 
may have to consider a disturbance of the ordered state 
by muons. 

We mention that Paixao et al. have proposed for the 
ordered state of Np02 that a triple-g r 5u -type octupole 
ordering is realized and a triplc-q r 5g quadrupole mo- 
ment is induced. 2 -* From the observed muon spin preces- 
sion frequency of to = 7 MHz 38 ) and the muon gyromag- 
netic ratio 7 P = 2ir x 13.55 MHz/kG at 8 K, we estimate 
the internal mag netic field as H = uj/j^ ~ 500 G. 39 ' 
This internal field is smaller than that of antifcrromag- 
netic UO2 by an order of magnitude. Thus the observed 
field in Np02 may be interpreted as coming from oc- 
tupole moments. However, the internal field in Np02 is 
still much larger than that of eq. (6.47). Provided the 
muon is really probing the intrinsic internal field of the 
systems, the muon stopping site should be about 1 A 
from a Np ion. 

7. Pressure Induced Moment 

It has been discussed that antifcrromagnctism should 
be induced in an antiferro-octupolar phase by uniaxial 
stress. 5 ) We apply this idea to the case of the T 5u order. 
The uniaxial pressure p along the [001] direction accom- 
panies the r 3 strain 



r. . (7.65) 

V3(Cu - C 12 )/2 

We solve the mean field equation with the finite strain 
e u and the corresponding quadrupole-strain interaction 
<7r 3 (sec, cq. (4.18)). The antifcrromagnetic moment is 
induced in the a;y-plane since (T^ u ) ^ together with 
(J 2 - J 2 ) 7^ gives (J x ) ^ 0. The direction of the anti- 
ferromagnetic moment is along [110] or [110]. In Figs. 9 
and 10, we show the pressure dependence of the octupole 
and magnetic moments at absolute zero with experimen- 
tal values: (Cn - C 12 )/2 ~ 20.4 x 10 11 erg/cm 3 and 
|.or 3 | = 120 K. 11 ) The mean field solution for (T z 5u ) be- 




1 1.5 2 2.5 
p(GPa) 



Fig. 9. Pressure dependence of the antiferro-octupole moment 
(T|«> (= (T5«» and (T|">. 



u 



V 




Fig. 10. Pressure dependence of the antifcrromagnetic moment 
Mf (= Mf). M Z AF = 0. 



comes zero for p > 0.9 GPa. This implies that both the 
octupole and dipole moments lie in the xy-planc. With 
p = 1 GPa, the magnetic moment is estimated to be 
0.88y(iB (0.82/ib), when gr 3 is positive (negative). This 
value should actually be reduced by quantum fluctua- 
tions neglected here. However, it is likely that this value 
remains within the experimental range. 

8. Summary 

To summarize, we have proposed that the T 5u octupole 
moment is a plausible candidate for the order parameter 
of the phase IV in CezLai-^Bg. This octupole ordered 
state is consistent with the lattice distortion along the 
[111] direction, 13 ' 14 - ) the softening of the elastic constant 
C44, 11 ) the broken time reversal symmetry, 16 ' 17 ) the ab- 
sence of dipole moment, 18 ^ and the cusp in the magne- 
tization at Ti.iv- 12 - 1 Especially, the sudden softening of 
C44 and the cusp in the magnetization at Tj.iy are diffi- 
cult to be explained by a pure quadrupole order, but are 
reproduced by the present octupole ordering model. 

Many properties of the phase IV are consistent with 
the octupole order, but the internal field estimated by our 
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mean field theory is much smaller than that suggested by 
the /xSR experiment. 17 ) We recall that the NMR spectra 
become broad in the phase IV, 16 ) and the spectra cannot 
be explained by the static and staggered octupole mo- 
ments, either. Hence some dynamical or quantum fluctu- 
ation effects may be playing a significant role. Clarifying 
dynamical aspects and identifying the periodicity in the 
phase IV are the most challenging open problems to be 
addressed in the near future. 

We have also shown that antiferromagnetic moment 
lying in the xy-plane is induced under uniaxial pressure 
along the z direction. The estimated value of the antifer- 
romagnetic moment is of the order of 1/iB: and we wait 
for experimental efforts to detect the moment. 
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